Electron Cotunneling in a Semiconductor Quantum Dot 
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We report transport measurements on a semiconductor 
quantum dot with a small number of confined electrons. In 
the Coulomb blockade regime, conduction is dominated by co- 
tunneling processes. These can be either elastic or inelastic, 
depending on whether they leave the dot in its ground state 
or drive it into an excited state, respectively. We are able to 
discriminate between these two contributions and show that 
inelastic events can occur only if the applied bias exceeds the 
lowest excitation energy. Implications to energy-level spec- 
troscopy are discussed. 

PACS numbers: 73.23.Hk, 73.40. Gk, 73.61.Ey 



Quantum-dot devices consist of a small electronic is- 
land connected by tunnel barriers to source and drain 
electrodes |Q. Due to on-site Coulomb repulsion, the 
addition of an electron to the island implies an energy 
change U — e 2 /C, where C is the total capacitance of 
the island. Hence the number of confined electrons is a 
well-defined integer, N, that can be controlled by vary- 
ing the voltage on a nearby gate electrode. Transport of 
electrons through the dot is allowed only at the transition 
points where the N- and (N + l)-states are both energet- 
ically accessible. Otherwise, N is constant and current is 
strongly suppressed. This is known as Coulomb blockade 
0. At low temperature, however, higher-order tunnel- 
ing events can become dominant. These arc commonly 
known as cotunneling events since they involve the simul- 
taneous tunneling of two or more electrons 0. Electron 
cotunneling has received considerable attention over the 
last decade. Initially it was recognized as a limitation to 
the accuracy of single-electron devices. More recently, it 
has acquired a broader relevance, especially since an in- 
creasing activity has been focused on quantum dots with 
a deliberately enhanced tunnel coupling to their leads. 
These systems allow the investigation of high-order trans- 
port processes and many-body phenomena, such as the 
Kondo effect H|}. The latter can be regarded as the 
result of a coherent superposition of higher-order cotun- 
neling events. Here, we will focus on the lowest order 
contribution to cotunneling. 

Previous experiments were performed with metallic is- 
lands or large semiconductor dots |^-|Tc|] , where the 
energy spectrum is essentially continuous and many lev- 
els contribute to cotunneling. Here, we study cotunneling 
through a small quantum dot where the energy levels are 
well separated, and where the absolute value of N is pre- 



cisely known. A cotunneling event is called inelastic when 
it leaves the dot in an excited state. Otherwise it is clas- 
sified as elastic. We identify two regimes: one consisting 
of elastic processes only, and one including both elastic 
and inelastic contributions. We note that the transition 
between these regimes can be sharper than the charac- 
teristic life-time broadening of the dot states. In such a 
case, the onset of inelastic cotunneling can be exploited 
to measure the energy spectrum of a quantum dot with 
improved resolution. 

The stability diagram of a generic quantum dot can 
be obtained by plotting the differential conductance 
(dl/dVsd) as a function of bias, V s d, and gate volt- 
age, V g . Coulomb blockade occurs within the diamond- 
shaped regions in Fig. fa. The diamond size is pro- 
portional to the addition energy, defined as E a( id(N) = 
Hdot(N + 1) - Hdot{N), where n dot (N) is the electrochem- 
ical potential of an iV-electron dot. Inside the A^-electron 
diamond, u dot (N) < n Ll u R < u dot (N +1), with fi L , u R 
the Fermi energies of the leads. The diamond edges cor- 
respond to level alignment: u dot (N) — uj, or ur (see an- 
gled solid lines) . This alignment determines the onset for 
first-order tunneling via the ground state of the dot, lead- 
ing to a peak in dl / dV s d(V s d) ■ The onset for first-order 
tunneling via the first excited state occurs at a somewhat 
higher bias (see dot-dashed lines in Fig. la, and the cor- 
responding diagrams in Fig. Ib and le). These first-order 
processes have been exploited as a spectroscopic tool on 
the discrete energy spectrum of dots ]jj . 

Here, we are interested in second-order tunneling of 
charge, which becomes more apparent when the tunnel 
coupling between the dot and the leads is enhanced. We 
neglect contributions from spin that could give rise to 
the Kondo effect. Elastic cotunneling is the dominant 
off-resonance process at low bias. It gives rise to cur- 
rent inside the Coulomb diamond (light-grey region in 
Fig. la). The corresponding two-electron process (Fig. 
lc) transfers one electron from the left to the right lead, 
thereby leaving the dot in the ground state. 

For e\V sd \ > A (AT), where A (AT) is the lowest on- 
site excitation energy for a constant N , similar two- 
electron processes can occur which drive the dot into an 
excited state. For instance, an electron can leave the dot 
from the ground state to the lowest Fermi sea, while an- 
other electron from the highest Fermi sea tunnels into the 
excited state (see Fig. Id). Although this type of pro- 
cess is called inelastic H , the total electron energy is con- 
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served. The on-site excitation is created at the expense 
of the energy drop eV s d- To first approximation, the on- 
set of inelastic cotunneling yields a step in dl / dV s d(V s d) 
This step occurs when e|V s d| = A (TV), which is not 
or only weakly affected by V g (see also Fig. lc in Ref. 
IP^I ). As a result, inelastic cotunneling turns on along 
the vertical (dotted) lines in Fig. la. At the edge of the 
Coulomb diamond the condition for the onset of inelastic 
cotunneling connects to that for the onset of first-order 
tunneling via an excited state (dot-dashed lines). 

Our device has the external shape of a 0.5-/j,m-high 
pillar with a 0.6 x 0.45 /jm 2 rectangular base (inset to 
Fig. 2). It is fabricated from an undoped AlGaAs(7 
nm)/InGaAs(12 nm)/AlGaAs(7 nm) double barrier het- 
erostructure, sandwiched between n-doped GaAs source 
and drain electrodes The quantum dot is formed 
within the InGaAs layer. The lateral confinement po- 
tential is close to that of an ellipse |14|. Its strength 
is tuned by a negative voltage, V g , applied to a metal 
gate surrounding the pillar. A dc bias voltage, V s d, ap- 
plied between source and drain, drives current vertically 
through the pillar. In addition, we apply a small bias 
modulation with rms amplitude V ac = 3 /iV at 17.7 Hz 
for lock-in detection. Measurements are carried out in a 
dilution refrigerator with a base temperature of 15 mK. 
We find an effective electron temperature T e = 25 ± 5 
mK, due to residual electrical noise. 

Figure 2 shows d//dV^ in grey-scale versus (V s d,V g ) 
at 15 mK. Diamond-shaped regions of low conductivity 
(light grey) identify the Coulomb blockade regimes for 
TV = 1 to 4. The diamonds are delimited by dark nar- 
row lines (dI/dV s d ~ e 2 //i) corresponding to the onset 
of first-order tunneling. For TV — 1, as well as for TV = 3, 
sub-gap transport is dominated by elastic cotunneling 
with no evidence for inelastic cotunneling. The differ- 
ential conductance is uniformly low inside the Coulomb 
diamond. (Slight modulations are seen due to a weak 
charging effect in the GaAs pillar above the dot Jl5|.) 
This is different for N — 2, where the onset of inelastic 
cotunneling is clearly observed. As argued before, this 
onset follows (dotted) lines, nearly parallel to the V g axis 
p6[ . At the diamond edges they connect to (dot-dashed) 
lines where first-order tunneling via an excited state sets 
in. Similar considerations apply to N = A. 

The different behavior observed for TV = even and TV = 
odd stems from the fact that inelastic cotunneling occurs 
only if E a dd(N) > A (TV), as apparent from Fig. la. In 
the case of non-interacting electrons, E a d d (N) = U(N) 
for TV = odd and E add {N) = [/(TV) + A (TV) for TV = 
even, where U (TV) is the charging energy for TV electrons; 
energy levels are spin degenerate and consecutively filled 
with pairs of electrons. This is a reasonable picture if 
the level spacing exceeds the exchange interaction energy 
jL7| . In our small quantum dot the first three levels are 
indeed widely spaced as already discussed in Ref. Q . Fol- 
low TV, A (TV) exceeds not only the exchange energy but 
also U(N). This implies that for TV = odd, A(TV) lies 
outside the Coulomb diamond (i.e. A(TV) > E ad d(N)) 



and thus inelastic cotunneling is not observed. (Note that 
cotunneling between spin-degenerate states is an elastic 
process as initial and final state have the same energy.) 
For TV = even, A (TV) is always smaller than E ad d(N) and 
inelastic cotunneling can be observed. 

We now discuss the difference in life-time broadening 
between first- and higher-order tunneling. At the onset of 
first-order tunneling a certain level is aligned to one of the 
Fermi energies. In this case, an electron can escape from 
the dot, which leads to a finite life-time broadening of the 
observed resonance by an amount hT. Here, T = Tl+Tr, 
where Tl and Tr arc the tunneling rates through the left 
and the right barrier, respectively. (Note that these rates 
are independent of V s d, since our bias window (~meV) 
is much smaller than the height of the AlGaAs tunnel 
barriers (~50 meV). 

The onset of inelastic cotunneling is also characterized 
by a width. In the zero-temperature limit, this is deter- 
mined by the life-time broadening of the excited state. 
Two types of situations can occur. First, the excited 
state can be between \il and (see right inset to Fig. 
3) so that inelastic cotunneling can be followed by first- 
order tunneling. Such a decay event leads to a life-time 
broadening of at least TiTr m KT/2. Second, the ground 
and excited state are both well below [i^ and fiR, imply- 
ing that only higher-order tunneling is allowed (see right 
inset to Fig. 4). Now, decay from the excited state can 
only rely on cotunneling. For this higher-order perturba- 
tion, the corresponding rate, r co , is much smaller than T, 
leading to a reduced life-time broadening. To illustrate 
these arguments, we select different dI/dV sd — vs — V s d 
traces and analyse their shape in detail. 

Figure 3 shows two traces for TV = 2, taken at 15 mK 
for gate voltages at the horizontal lines in the left in- 
set. The dashed trace has several peaks. The two inner 
ones, at \V s d\ « 1.1 mV, correspond to first-order tunnel- 
ing of the 3rd electron via the 3-electron ground state; 
i.e. /irf t(3) = )1l or hr. The right (left) peak has a full 
width at half maximum (FWHM) of «200 («400) /iV. 
This is somewhat larger than the width, KT/e ~ 150 fxV, 
measured in the zero-bias limit. Indeed a finite bias al- 
lows non-energy-conserving tunneling events leading to 
additional broadening. The most likely source for energy 
relaxation is acoustic-phonon emission |l8j. The follow- 
ing pair of peaks, at \V sd \ « 2 mV, corresponds to the 
onset of first-order tunneling via the first excited state for 
TV = 2 (see Fig. lb). Because of the larger bias voltage, 
these peaks are visibly broader than the inner ones. Ad- 
ditional peak structures occur near the edges of the bias 
window. The origin of these peaks can not be precisely 
identified. 

The solid trace contains structure from both first- 
and second-order tunneling. The peaks labeled by open 
squares arise from first-order tunneling at the edges of the 
Coulomb diamond (see left inset). Steps, labeled by open 
circles, identify the onset of inelastic cotunneling and cor- 
respond to the open circles in the left inset. Their differ- 
ent heights are probably due to a left-right asymmetry 
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in the tunnel coupling to the leads. Their V^-position, 
which is symmetric around zero bias, provides a direct 
measure of A(2). The width of these steps is «150 fiY 
fl9|| . Since A(2) « U(2), the first excited state lies un- 
avoidably within the bias window when \V s d\ — A(2)/e 
and hence is allowed to decay into the lowest-energy lead 
(see the right inset to Fig. 3). As argued above, this sit- 
uation leads to a step-width exceeding hTn/e, consistent 
with our finding. Another structure occurs at V s d ~ 2.6 
mV and is probably due to the onset of inelastic cotun- 
neling via the second excited state for N = 2. The cor- 
responding line in the stability diagram is hardly visible 
due to its vicinity to the diamond edge. 

To study inelastic cotunneling when both ground and 
excited state lie well below the Fermi energies of the leads 
(Fig. Id) we need A(N) < E add (N). To this aim we 
move to N = 6, since A (6) can be effectively tuned by a 
magnetic field applied along the vertical axis. We tune 
the field to 0.35 T, such that A(6) w 0.1 meV, i.e. several 
times smaller than E a dd{6)- From a previous study we 
know that the ground state is a spin singlet, and the first 
excited state is a spin triplet pj|. The dI/dV s d — vs — V sd 
traces shown in Fig. 4 are taken at two different tempera- 
tures, but for the same V g , at the horizontal line in the left 
inset. The solid trace (15 mK) shows a broad minimum 
around V sd — 0, where transport is dominated by elastic 
cotunneling via the ground state (see also the light-grey 
region in the left inset). The differential conductance 
increases rapidly at the onset of inelastic cotunneling 
with a step- width of «20 fiY, i.e. much smaller TiT/e. 
This reduced width stems from the fact that the excited 
state can not decay directly into the lower energy lead 
(sec right inset). The corresponding life-time broadening, 
hT co , can be estimated from the cotunneling current, I co , 
at V sd = A(6)/e. We find ?ir co = M c0 /e « (ft/e) J A(6)/e 
dJ I dVsd(V s d) dV s d — 10 ^eV, consistent with the ob- 
served step-width. At T e = 25 mK the thermal broaden- 
ing of the Fermi distribution functions leads to a step- 
width of 5.Uk B T e /e ~ 12 ^eV |o|. Hence life-time 
broadening has been reduced here to the thermal limit. 

The cotunneling onset in Fig. 4 shows a peak struc- 
ture at low temperature (solid trace) in addition to the 
expected step structure. This is likely due to Kondo cor- 
relations, as discussed in Ref. On increasing temper- 
ature to 200 mK these Kondo correlations are suppressed 
such that only lowest order cotunneling contributes. This 
recovers the step structure (dashed trace). 
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FIG. 1. (a) Stability diagram in the plane of (V„d, V g ). 
Angled lines correspond to alignment of a dot-state with the 
Fermi energy of the leads. In this case, first-order tunneling 
sets in, or is increased, as illustrated in (b) and (e). In the 
light-grey area in (a), conduction is due to elastic cotunneling 
via virtual events as shown in (c). For eV s d > A(N), inelastic 
processes, illustrated in (d), increase the cotunneling current 
(dark-grey areas). A(N) is the energy spacing between the 
ground state and the first excited state, which in (b)-(e) are 
represented by solid and dotted lines, respectively. 
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FIG. 2. Measured stability diagram of our quantum dot 
at 15 mK and zero magnetic field. dI/dV s d is plotted in grey 
scale as a function of (V s d, V g ). Dotted lines have been super- 
imposed to highlight the onset of inelastic cotunneling. The 
dot-dashed lines indicate the onset of first-order tunneling via 
an excited state. Inset: scanning electron micrograph of the 
device. 

FIG. 3. Differential conductance as a function of bias for 
V g = -1.40 V (solid line) and V a = -1.30 V (dotted line) 
at 15 mK. These traces are extracted from the stability di- 
agram shown in the left inset. The horizontal lines indicate 
the corresponding V g values. The right inset shows the quali- 
tative energy diagram corresponding to the onset of inelastic 
cotunneling for N = 2. The horizontal arrow represents the 
possibility for an electron in the excited state to decay directly 
into the right lead by first-order tunneling. 

FIG. 4. Differential conductance as a function of bias at 
V g = —0.685 V. (Note that the bias window is much smaller 
than in Fig. 3.) The solid (dashed) line is taken at 15 mK 
(200 mK). Left inset: stability diagram at 15 mK, around 
the 6-electron Coulomb diamond. The horizontal line is at 
V g — —0.685 V. Right inset: qualitative energy diagram cor- 
responding to the onset of inelastic cotunneling for N = 6. 
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